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Abstract 

We study Higgs diphoton decays, in which both photons undergo nuclear conversion to electron- 
positron pairs. The kinematic distribution of the two electron-positron pairs may be used to probe 
the CP violating (CPV) coupling of the Higgs to photons, that may be produced by new physics. 
Detecting CPV in this manner requires interference between the spin-polarized helicity amplitudes 
for both conversions. We derive leading order, analytic forms for these amplitudes. In turn, we 
obtain compact, leading-order expressions for the full process rate. While performing experiments 
involving photon conversions may be challenging, we use the results of our analysis to construct 
experimental cuts on certain observables that may enhance sensitivity to CPV. We show that there 
exist regions of phase space on which sensitivity to CPV is of order unity. The statistical sensitivity 


of these cuts are verified numerically, using dedicated Monte-Carlo simulations. 
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I. INTRODUCTION 


The recent discovery of a Higgs-like boson [1, 2] has prompted intense interest in the 
precise measurement of its couplings and properties. Such measurements are a direct probe 
of new physics (NP) beyond the Standard Model (SM), especially since many extensions of 
the SM modify the Higgs couplings to gauge bosons and fermions. 

Of particular interest is a search for parity and CP violating Higgs decays, since these 
would be a clear signal of NP [3-20]. In Higgs decays to vector bosons the CP violating 
effects can only be due to irrelevant NP operators. Fortunately, in the SM h > yy (and 
h — Zy) decays are also due to irrelevant operators, with the first non-zero contribution 
occurring at one-loop. Thus we can expect large CP violating effects from weak scale NP 
inh + yy. In contrast, the h > ZZ* and h 5 WW* decays proceed in the SM through 
relevant tree-level operators tightly related to the Z and W masses. CP violating effects 
from NP are expected to be comparatively small in these decay modes. 

In the presence of CPV, the total h — yy decay rate must be proportional to the sum of 
squares of CP-even and CP-odd terms — i.e |CPeven|?-+|CPoaa|? — and therefore, by comparing 
the h — yy rate to the SM expectation, one may probe for NP directly. However, this type of 
search cannot distinguish between CP-even and CP-odd NP contributions to the total rate. 
Moreover if the CP-odd contribution is small, then CPV signals are quadratically suppressed, 
and if it so happens that NP enters both the CP even and odd terms, such that the total 
h — yy rate matches the SM expectation, then one cannot detect NP at all. Probing 
the differential h — yy rate for CPV ameliorates these problems. In the first instance, 
the differential rate may feature an interference term of the form 2CPeyenCPoag. Combined 


with non-interfering terms, one may distinguish CP-even and CP-odd NP contributions. 
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FIG. 1. An illustration of an example of a CPV sensitive observable in h + yy — 4e. The 
Higgs decays to on-shell photons which convert in the detector. The distribution of the azimuthal 
angle y between the two planes formed by each positron and its parent photon depends on the 
Higgs couplings to CP even and odd operators. The electrons do not need to be co-planar with 
the corresponding photon-positron planes. The positron-photon plane is shown in magenta and 


the electron-photon plane in blue. For further details and subtleties see the main text. 


Secondly, small CPV signals are only linearly suppressed in this interference term. 


The h — yy phase space distribution alone, however, is not sensitive to CP violating 
effects, since the Higgs decays isotropically to two photons. Nevertheless, the underlying 
CPV structure in the differential h — yy rate may be determined if one is able to measure 
the linear polarizations of the outgoing photons. This in itself is an old idea, first proposed 
for the determination of the 7° parity [21-23] and, more recently, to probe NP effects in 
radiative B decays [24]. It relies on the fact that a spin-0 particle decays to either two positive 
or two negative helicity photon states, which acquire a relative CPV phase in the presence 
of non-trivial CP structure. The linearly polarized photon states are a superposition of both 
helicities, permitting one to extract this CPV phase. It is not feasible to directly measure the 
linear polarization of O(60 GeV) photons from Higgs decay. However, in both the ATLAS 
[25] and CMS [26] detectors roughly half of the photons from Higgs decays convert via the 
well-known Bethe-Heitler (BH) process into ete” pairs inside the silicon tracker. This has 
an important benefit: the orientation of the produced ete” pairs encodes the underlying 


CP structure. Figure 1 illustrates an observable expected to be sensitive to CPV. 


Previous proposals to measure CPV in h — yy, or in other neutral meson diphoton 


decays, via double photon conversion appear in Refs. (3, 23, 27]. We extend these studies 
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by performing an analysis of the actual manner in which the leptonic phase space encodes 
the CP violating effects. We examine the encoding of CPV in the doubly converted Higgs 
diphoton decay — hereafter the Higgs-Bethe-Heitler (HBH) process — both analytically and 
numerically, making use of the spinor-helicity formalism to obtain compact expressions for 
the full differential scattering rate and its leading order terms. 

A key difficulty in extracting CPV signals from differential scattering rates like dT'/dy, 
defined in Fig. 1, is that the signal is largely washed out under integration over the other 
phase space variables. However, in this work we use our analytic control of the full dif- 
ferential scattering rate to show that large CPV signals may be achieved. We find certain 
observables exhibit O(1) sensitivity to CP violating effects on small regions of the phase 
space, corresponding to a small fraction of the converted decays. These regions are identi- 
fied by sensitivity parameters, derived from our analysis. We show analytically and confirm 
numerically that if one cuts the HBH event data according to these sensitivity parameters, 
the CPV signal is dramatically improved. 

Performing such an experiment will be difficult. For example, one must resolve the 
opening angles in the photon conversion, and this requires extremely accurate tracking 
resolution. Other subtle effects, such as soft scatterings of the electron and positron in 
the detector material need also to be examined. In addition, reaching a sensitivity to Higgs 
couplings would require obtaining a signal-rich sample of h + yy events. Our approach here 
is to defer these considerations to future work, and consider mainly the theoretical aspects of 
this process. In doing this, our intent is to motivate a very challenging measurement, perhaps 
to be done after the LHC upgrade or in a future Higgs factory. We note however, that as 
a warm-up to Higgs studies, photon conversions can be studied in background samples as a 
test of the standard model (we present the phase space structure of gg — yy with converted 
photons in an appendix). 

This paper is structured as follows. In Section II we first present a motivating phe- 
nomenological analysis of the expected size of CP violation in Higgs decays, followed by an 
analysis of the unconverted h — yy process. This is presented in the language of helic- 
ity amplitudes and we show that the CPV terms in the differential rate arise from helicity 
amplitude interference. In Section III we discuss photon propagation and conversion in mod- 
ern detectors, angular resolution limits, and the central role of the nuclear form factor in 


Bethe-Heitler conversion. In Section IV we then proceed to examine the HBH process itself, 


presenting explicit results within the spinor-helicity formalism for the leading order HBH 
square amplitude. This is followed in Section V by a derivation of CPV observables and 
their sensitivity parameters for the special case that the Higgs is at rest. These sensitivity 
parameters can be used as cuts, which extract the phase space regions on which we expect 
O(1) CPV effects. Numerical simulations are presented which confirm these expectations, 
and further compare the performance of the various sensitivity parameter cuts. In this work 
we focus on the h — yy decay, but our analysis can be used to examine other searches, e.g., 


for h + Zy or other decays involving converted photons. 


II. HIGGS DIPHOTON DECAY WITH CP VIOLATION 
A. Motivation for measuring CP-violation in di-photons 


The CP structure of Higgs decays was already studied experimentally in h + ZZ” decays 
(28, 29], with pure pseudo-scalar coupling disfavored at more than 30 [29]. Still, there is 
strong motivation to measure CPV also in loop induced h — yy decays. The motivation 
is based on the expected sizes of CP violating and CP conserving terms in h > yy and 
h => ZZ* decays. After electroweak symmetry breaking (EWSB), the relevant terms in the 


effective Lagrangian at the scale u ~ m,/2 are 


[eee ENI, +c REF” Fezz- hZ" Zyy 
v TU TU 

+ zzz ħZuŽ" + Poia (1) 
in which Fr, and Z,, are respectively the photon and Z field strengths, and XH = rrop Xag, 
e123 — 1, while v = 246 GeV is the Higgs vev. Taking Higgs to be a scalar, the first line 
of Eq. (1) contains CP even and the second line CP odd operators.! Present data imply 
cy = 1.04 0.13 [30], assuming CP conservation. If either čzz or č are found to be non- 
zero,” CPV in Higgs couplings and thus NP will be discovered. The couplings c and č arise 


at one-loop in perturbative UV theories, and can be at most O(1) in order to agree with 


the observed h — yy rate. For example, č can arise from a massive NP fermion loop that 


1 Since all the couplings are C even, P and CP violation are the same. While we hereafter always refer to 


CP violation, it should be kept in mind that it is equivalent to parity violation for this effective theory. 
2 In fact, č may be non-zero at three loops in the SM [31]. However, since this contribution falls well 


below the feasible detection threshold for the experiments under consideration in this work, we neglect 


this contribution, and treat the SM value as č = 0. 


T 


is axially coupled to the Higgs. Note that in Eq. (1) we integrated out both W and t loop 
contributions to c, so that csm = —0.81 in the SM. Similarly, the dimension 5 couplings 
of Higgs to ZZ, czz and čzz, also arise at one-loop order. In generic NP models we thus 
expect c ~ czz and Č ~ Gzz. 

The h + ZZ* decay is dominated by the CP even renormalizable coupling cy, while the 
h — yy decay is given by higher dimensional operators. The relative size of CP violating 
effects in any channel is given by the ratio of the interference terms to the total amplitude 
squared. For h + ZZ* the CP odd interference term is proportional to (a/27)ézzcy, while 
the total squared amplitude is dominated by c&. The typical size of CPV observables in 
h — ZZ* is therefore set by the ratio of the two, 


a Ba, 3 O(10~*), (2) 


for O(1) couplings. In the diphoton channel both terms are loop suppressed and the figure 


of merit for CP violation is i 
T5 aa ~ 01), (3) 
again assuming O(1) couplings. As will become clear, the measurement of CPV in h > yy 
is a challenging one, especially in comparison to the relatively straightforward measurement 
in h + ZZ*. However, the expected size of the effect may partially compensate for this. 
In addition to the CPV observables discussed in this paper, the CPV operator in Eq. (1) 
also modifies the overall h — yy decay rate, so that 
. Thy) e+ 
nT hys Bm 


i (4) 


where [(h — yy)sm is the SM rate. The total rate is only quadratically sensitive to CP 
violating NP because the interference terms integrate to zero over phase space.In contrast, 
the differential rates contain CP odd terms proportional to cc and thus may be linear in the 
NP coupling. This can lead to substantial increase in sensitivity for small č. 

Before proceeding, two remarks are in order. First, it is important to mention that there 
are severe constraints on ye€ from bounds on the electric dipole moment of the electron [32— 
34]. Taking the electron yukawa, ye, to be the SM one, this gives č < 107°. These bounds are, 
however, absent in the limit where the 125 GeV Higgs does not couple to electrons (for other 
possibilities in concrete UV models, see (32, 35]). A strong motivation for contemplating a 


non-zero value of ¢ is, for instance, that it would be generated by new CP sources in models 


that lead to electroweak baryogenesis, see e.g., [35]. An independent measurement of č is 
thus desirable. 

Second, we assumed here that c and č are real. In the SM c obtains its dominant 
contribution from a W loop and a smaller destructive contribution from a top quark loop. 
However, there is also a smaller contribution from b quark and light quarks. These can go 
on-shell in the loop, generating complex effective couplings c and č. This means that c and 
č obtain a relative strong phase. The result of such a strong phase, when combined with the 
weak phase, would be to induce direct CPV in decay, such that the decay rate of the Higgs 
into two positive helicity photons is not the same as into two negative ones. These strong 
phases are of ~ O(1%) in the SM (36, 37], and we assume the strong phases are similarly 
small for NP effects. Consequently, we neglect direct CPV, and assume that c and č are 


real. 


B. Helicity interference 


We proceed to examine the unconverted h — yy process. The effective operator mediat- 
ing h — yy decay has the general form, cf. Eq. (1), 


a Ca ~ 
hF F” — -— hF F" . 5 
y ” 2nv ” (5) 


T T 


Heg = —cC 


For a Higgs that is a scalar, the first term is CP even and the second is CP odd. CP is 
therefore violated if the CP phase 


Esta VGA , (6) 


is found to be non-zero. 


The h — yy helicity amplitudes are (dropping the overall a/mv factor, cf. Eq. (5)) 


ka, AM 


= | (Ry KEY (SY) — ha (Ea ey Aka, (E), ka CRY] , 
(7) 
where k; are the photon momenta, £ is the polarization vector of the ith photon (i = 1,2) 


with helicity A; = +, and e[p,q,r,s] = puqrpsceh’’’. A Latin subscript i = 1,2 hereafter 
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denotes the corresponding photon. To compute helicity amplitudes we employ the spinor- 
helicity formalism, see Appendix A for our conventions and a brief review. Using Eqs. (A5), 


(A6) and (A10), one finds that the non-zero helicity amplitudes are 


ME = m2 (c+ id) = mV 24+ Cet, (8) 
while M=* = 0 as expected from angular momentum conservation (cf. also the results of, 


e.g., (38, 39] for h + ZZ). We see that CPV introduces a relative phase, €, between the 
two-photon helicity amplitudes. Furthermore, a differential rate may depend on € if and 
only if there is interference between Mtt and M77, or more precisely, between amplitudes 
involving the ++ and —— photon helicity configurations. We call such interference helicity 
interference. 

Let us now translate Eq. (8) into a Hilbert space language. The final states of h > yy 
decay are the two-photon states |+-+) and |— —), with + indicating the helicity of each 
photon, so that CP|4+) = |F F). Eq. (8) then translates to 


Healh) x VR +(e |++) +e ]--)) (9) 


As above, the CP phase € appears as a relative phase between the |++) and |——) terms. 


Now, the total rate for h — yy decay is proportional to 


5 |(flHerlh)P = So MI. (10) 


fat+-— Tari 


Orthogonality of |++) and |— —) ensures that the total rate is independent of €, i.e., there 
is no helicity interference. In contrast, any experiment for which the final state is a linear 
superposition of the two helicity states would generate helicity interference. This is the case 
at collider experiments in which the on-shell photons with definite helicity are intermediate 
states: The final state is a converted photon — an ete” pair with a particular set of momenta 


— which has non-vanishing overlap with both helicities. 


C. A thought experiment with polarizers 


The overlap of each photon helicity with a BH pair will determine the strength of helicity 


interference and our ultimate sensitivity to the Higgs CP properties. The details at the 
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FIG. 2. A linear polarization thought experiment in Higgs rest frame. Pi are linear polarizers 
oriented orthogonal to the photon momentum direction. The angle ¢ is measured between the 
linear polarization vectors e's. 

level we need are quite involved, and will be described in Sec. IV. As a warm-up we instead 
consider a thought experiment in which we can measure linearly polarized photons. 

Let us imagine that we have been able to manufacture a linear polarizer for gamma 
rays. We produce a Higgs at rest between polarizers Pı and P5», such that each photon 
travels through a polarizer (see Fig. 2) before being absorbed by a detector that counts 
photons. The polarizer P; (i = 1,2) projects an incoming photon onto a linearly polarized 


state, |d;) = e*|+) + ehi 


—), that has polarization oriented at angle ¢;. From Eq. (9) the 


amplitude of the two-photon wave function observed by the detectors is 
(eta) @ (H+ (1) Han) ~ Ve +P colo + £) , (11) 


where 6 = &1 — His a relative azimuthal angle between the two polarizers. As the angle ¢ 
is changed, the differential rate in the detectors changes as cos?(¢ + £). One finds 

dl’ 2 a? m? 7 

dd = zon cos (P4), Dhow = Tag PE) (12) 
Note that the CP odd term in the differential rate (12) is proportional to sin 2€ sin2¢ œ 
cčsin 2ġ. The differential rate is thus linearly sensitive to CPV coupling č, whereas the total 
rate is quadratically sensitive. 

In summary, we have shown that only the terms receiving contributions from both of the 


definite helicity two-photon amplitudes, so that there is helicity interference, are sensitive 


to the CPV phase €. These interference effects can in principle be of O(1) in size. 
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FIG. 3. The contributions to photon cross-section on 78Si, Ty (?8Si), from BH ete- pair production 
in nuclear field (solid blue line), pair production due to scattering on electron cloud (red dashed), 
Compton scattering (dot-dashed yellow) and Rayleigh scattering (magenta double dot-dashed), as 


a function of photon energy Ey. Calculated using NIST’s XCOM database [40]. 
III. BETHE-HEITLER PHOTON CONVERSION 


We now study the process that can be used for photon polarization measurement, namely 
the conversion of a photon into an ete” pair in matter. In this section we study a conversion 
of single isolated photon, which we will then use in the Sec. IV for the case of h > yy with 


converted photons. 


A. Photon propagation and conversion 


Photon conversion to ete~ pairs may proceed either by Dalitz conversion in vacuum, 
for an off-shell intermediate photon, or by BH conversion on atomic nuclei, which occurs 
for on-shell photons (for a review see, e.g., (41, 42]). The Dalitz conversion rate carries a 
suppression factor of O(10~*), and is not of immediate practical interest. Moreover it mainly 
proceeds via a longitudinal photon so that the above helicity analysis no longer applies. In 
contrast, the CMS and ATLAS pixel detectors contain a significant amount of material, so 
that ~ 50% of photons convert inside the tracking systems via the BH process [43, 44]. 
Based on the composition of the detectors, we assume in this work that the target nucleus 


is always 2Si, at rest in the laboratory frame. This nucleus is spin-0, and has no nuclear 
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FIG. 4. Cumulative distributions of Bethe-Heitler conversion cross-section for a photon with 
E, = 60 [GeV] scattering on a 28Si nucleus, with respect to various opening angle cuts. Three 


distributions are shown: P(Oge > Ocut), i.e. with photon-lepton opening angles, 6, unconstrained 


(blue line); P(A, and ge > Ocut) with electron-photon opening angle 6_ unconstrained (red dashed 


line): P(0x and ee > cut) (black dot-dashed line). 


magnetic moment. We therefore do not consider the effects of target polarization on the BH 
process [45]. 

One might be concerned by the prospect of photon polarization decoherence for the 
photons propagating inside the silicon. However, at photon energies ~ O(m,/2) the pair 
production in the nuclear field is by far the largest contribution to the photon scattering 
cross-section in an atomic lattice [46], see Fig. 3. As a result, to an excellent approximation, 
the photons remain coherent up until their BH conversion. We shall also assume that the 
BH scattering is quasi-elastic, i.e. that the target nucleus remains in a coherent state during 
and after the scattering. For the kinematics considered, the quasi-elastic limit is an excellent 


approximation of the full BH conversion [41, 42]. 


B. Angular resolution limitations 


Following the h — yy experiment discussed in Sec. IIC, in order to measure the CPV 
in a doubly converted h — yy decay, we might expect that angular distributions between 


planes formed by spatial momenta need to be measured. There are several possible planar 
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distributions that can be constructed, involving either: 


1. The ete~ plane formed by the ete momenta, or; 


2. The ye* plane formed by a lepton and its parent photon, as in Fig. 1. 


The first requires resolving the orientation of the leptonic spatial momenta. The second 
requires the orientation of the leptons with respect to their parent photon, which could be 
achieved by identifying the vertex associated with the Higgs (from other tracks in the event) 
as well as the location of the photon conversion, giving the photon direction. 

Such measurements require exquisitely precise tracking. Because the momentum transfer 
to the nucleus is small, the relative angular orientations between the photon and leptons are 
tiny in typical photon conversion: for mass m and energy E, the angular scale is typically 
m/E ~ 10“ for a 60 GeV photon conversion to electrons. There is however a distribution 
for these angles. In the limit of very large statistics one can hope to get a sample of events 
where these angles can be measured. 

In the ATLAS detector, for instance, the intrinsic accuracy in silicon pixels located be- 
tween 5cm and 12 cm from the interaction point is 10um in R — 6 direction and 115um in 
z direction. The intrinsic accuracy of SCT strips located between 30cm to 51cm from the 
interactions point is 17um (R — $) and 580um(z) [47]. One may therefore hope to measure 
the orientations of the ete” plane even for opening angles as small as 6e ~ 10” — 107°, 
where the relative leptonic angle 0, is defined by 
P+: P- 


i a (13) 
Ip: Ilp-| 


cos Ope = 


for leptonic spatial momenta p+. By comparison, a 60 GeV photon converting to a eet 


pair has an opening angle of Ore > 1074 in 38% of the cases and O > 10? in 4% of the 
cases. The full cumulative distribution for Oee > Oeut, is Shown in Fig. 4. 
Another experimental challenge is the multiple scattering of outgoing electrons when 


traversing the detector medium. This can affect the measurement of the electron direction 


and thus the orientation of the ete~ or ye“ plane. Using Eq. (30.15) in [46], the width 
of the angular distribution is ~ 1074 for a 30 GeV electron, assuming it traverses 0 = 0.1 
radiation lengths of the material. This width roughly scales as V¢/E, where E is the lepton 
energy. The measurement of polarization planes in the current and future detectors will 


thus be challenging, but may be achievable on a statistical basis. Bearing in mind these 
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FIG. 5. The elastic form factor G§!(q?). The dashed lines show the limiting behavior G5’ ~ a‘q* 
for |a? q?| < 1 (green dashed line) and G$ ~ 1 for |a? q?| >> 1 (red dashed line). The scale at which 
screening of the nucleus becomes important is denoted by a, which is smaller than the Si atomic 
radius, Tatom: At scales well outside the atom, corresponding to small —q?, nuclear conversion is 


suppressed by the form factor screening. 


experimental questions, in the remainder of this paper we adopt a theoretical approach to 
this problem: We consider a thought experiment where all angles can be resolved and explore 


the sensitivity to Higgs parameters in this best-case scenario. 


C. Nuclear form factor 


The BH conversion depends on a momentum transfer, q”, between the photon and the 
nucleus. Assuming quasi-elastic scattering, the photo-nuclear scattering is encoded in an 
elastic nuclear form factor G$![q?] (see Eq. (22) below). This form factor plays an important 
role of suppressing scattering at low —q?, that is, at scales larger than the Si atom. 

Let us discuss briefly the behavior of the nuclear form factor. The threshold for an 
E, ~ 60 GeV photon to convert to an ete" pair is at —q? = 4m* | B2 ~ 10718 GeV?, with 
m the electronic mass, but can occur at momentum transfers as large as —g? ~ 108 GeV?. 
This should be compared with the radius of the Si atom, ratom œ= 1.1A [48] or 72, œ 107? 
GeV’, and with the nuclear radius rpucear ~ 4 fm which gives ry2, ~ 1073 GeV?. Within 


the —q? range of interest for conversion — 10~'8 up to 1076 GeV? — the nuclear charge is thus 
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screened at low —q? by the atomic shell electrons. In this work we use the simple expression 
for the atomic form factor [49] described in detail in [41] and given by 
2,44 

CH) = banya (14) 
where a = 184.15(2.718)~'/?Z-/3/m and Z is the atomic number of the nucleus with mass 
M. For Si, a~? = 1.22 x 1019 GeV’. There are two limits of interest. The first is 
la? q?| >> 1 in which GS ~ 1, the second is the limit |a? q?| < 1 in which Gs! ~ qf, see Fig. 
5. That is, the form factor suppresses the BH cross-section for small —q? < a?. To the 
extent that the 1/g' factor in the BH cross-section determines the dominant phase space 
configurations of the final states, this suppression significantly alters the important regions 


of phase space for BH conversion up to the a~? scale. Specifically, the form factor increases 


the probability of significantly acoplanar BH conversions. 


IV. THE HIGGS-BETHE-HEITLER PROCESS 


In this section we present a formal analysis of the Higgs-Bethe-Heitler (HBH) process, 
h — yy with both photons converting to ete” pairs. The main result is that we obtain 
compact, leading order expressions for the HBH rate, and gain insight into the structure of 


the terms sensitive to CP violation. 


A. Amplitude and cross-section 


The amplitude of interest is given by a menorah diagram, consisting of a h > yy and 


two BH conversion subdiagrams, summed over internal photon polarizations, viz. 


P! 
ji 1 
P 
Nan 
P1.,71 
Pı}, S1 
(MARY, = Sarees kr Ar + lepton exchanges 
rs2rs = ARN. ko Ap i 
P24, 52 
p2_,72 
q2 
5 vA 
v 
P, 


= i D> [ele kog — kk) + (tro, P) (E (EF) [MEn], Mka], » 


A1,A2 


(15) 


in which A = 7vu/e®a, for QED coupling e. The BH amplitudes are 


Ple), a 


a 
yf yl 
: — Mm 


[Mul], = a” (p:_) A 


1 
Ki — P; 5, — 
where we have not yet taken nuclear form factors into account, and kept explicit the cor- 
responding Lorentz index of the nuclear electromagnetic current. The Latin subscripts, 
i = 1,2, label each photon, while s;, ri = 1,2 are respectively the positron and electron 
spins (see Eq. (A8)), and Aj = + the outgoing photon helicities from the h — yy vertex. 
We have also suppressed repetitions of the photon index, such that X;,, is henceforth a 


shorthand for X; We shall often refer to the BH subdiagrams for each photon as the 


photon branch. 

We assume both nuclei are initially at rest in the lab frame, so P” = (M,0) where M 
is the mass of the nucleus. The Higgs need not be at rest in the lab frame. As discussed 
in Sec. IITA, we assume quasi-elastic scattering, that is P? = P? = M?. This implies that 


2q:: P +q? = 0, and so 
E,-E;, —E, #@/2M —0, (17) 


where F;, Ei, are the photon and lepton energies respectively. It is also convenient to define 


1 qi: P 1 qi 
K5 pe _ gl = — | P” 4+ 18 
! al Ti 4 ) al t5 i ( ) 


under guasi-scattering conditions. The exchange energy with the nucleus q? = P? — P? « 


M. Le. the nucleus velocity is non-relativistic, so to an excellent approximation, it follows 


that 
Bye E, +E, Ge, (19) 


in the lab frame. 
We define the BH nuclear form factor tensors on each photon branch [41], 


p 
wee = —W,(q?) G = a + Wa (p)@0: , (20) 


1 
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such that the unpolarized HBH squared amplitude 


MP || 206m? — Pp?) = 2M E o MÄE Win Wav - (21) 

i=1,2 TiSi 
On the left of Eq. (21) we have factored out the d—functions that enforce quasi-elastic 
scattering. The form factor Wı(q?) = 0 for quasi-elastic scattering on a spin-0 nuclear 


target [45], while 
Wa(¢) = 2Mô(M? — P”)G3(q’) , (22) 


in which G$!(q’) is the form factor given in Eq. (14). Hence the unpolarized HBH squared 


amplitude reduces to 


s (q2) Meee 
MP = ng S N aM e Oud Q do . (23) 


TiSi 
B. Helicity structure 


One may calculate |M|? via the usual Feynman method, which relies on polarization 
completeness relations to compute traces. This approach leads to thousands of terms, of 
which the naively dominant terms cancel due to Ward identities. Extracting leading order 
expressions is therefore difficult, and moreover, high numerical precision is required for 
numerical stability. As an alternative, we employ a spin and helicity analysis combined with 
the spinor-helicity formalism to compute the HBH amplitudes. These may be subsequently 
squared and summed over external spins to produce the full HBH rate. In the following we 
provide a brief overview of these results, while details are provided in Appendices A and B. 


The BH spin-helicity amplitudes are defined as 


ajs, = [Minli Qin/G - (24) 
With reference to Eqs. (15) and (23), the HBH spin-helicity amplitudes are correspondingly 


VEE. (iy «gh? — KERS) tks) AN naa 3 


ea aa 
Tsers 
where we introduced the abbreviation 


G; = G (a) - (26) 


18 


Making use of the spinor-helicity formalism (see Appendix A), we obtain the spin amplitudes 


m? otf ds gs I 
Mirades = >` Mi? , = ay GiGoVe + €? (Pokan +e Kari.) i (27) 


A1à2 


These are reminiscent of equation (9) with final leptonic spin states (r;s;| and a; œ (s;r;|+), 


as expected. The HBH square-amplitude is correspondingly 


MP = » do Mug, = mo t= as > lot of. eS tar rs ra me f (28) 


1rs2rs NM A2 1rs2rs 


2 


Note that the sum over photon polarizations is inside the absolute value as expected for 


entangled h — 2y — 2(ete”) decay, but the sum over lepton spins is incoherent. Eq. (28) 


shows that the BH spin-helicity amplitudes aj” are all one needs to determine the entire HBH 
square amplitude. Parity invariance of the BH amplitudes relates amplitudes of opposite 


helicity and spins to their complex conjugates, 
on = NrsQrs » (29) 


where 5 is the opposite spin to s, and 72, = 1, Ms = Nrs- Hence we need only determine 


aj. We shall see below that spinor-helicity methods, when applied to af, also allow for a 
well-controlled expansion of dominant, sub-dominant and negligible terms in |M|?. 
We see explicitly from Eq. (28) that the helicity interference terms are equivalent to the 


£-dependent terms, as discussed in sections IIB and IIC. For a particular leptonic spin 


+ 


configuration {r;, si}, helicity interference occurs so long as aj and ay are both non-zero. 


However, we see in appendix B, and in particular in Eq. (B3), that for our particular choice 


of spinor basis (see Eq. (A8)) 


lanl = lai) ~ lar, > ls, | 7 (30) 


or equivalently Jas, | ~ laf,| ~ la, | lah, |. This hierarchy means that we may therefore 


discard any terms containing either af, or az as subleading. It follows that the leading 


122 111 
2 
+2 2 =i] 
ri Asi 


z >. ay, a3, e $+ A, Wana E 
nine ra #82 
(31) 


order squared amplitude is 


(CLC 


IM|? ~ m Mp AA2 Haa laman 


Ain Akr 
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while the leading interference term is 


c + Z i 
AA a5, Ay A5, i s (32) 


IMP. 2m, 


We see that only ay» (or equivalently 091 12) amplitudes enter the leading order £- 
dependent interference terms. 


The CP odd helicity interference term in (32) is proportional to cc, 


MP litical = a A T : (33) 
Note that interference terms between amplitudes, produced by the FF and FF operators 
respectively, are CP odd. However, the helicity amplitudes under consideration here receive 
contributions from both CP odd and CP even operators — manifestly they depend on € — and 
hence helicity interference terms contain both CP-even and CP-odd pieces. Consequently, 
we interpret the remaining piece of the helicity interference term to be the CP-even piece, 


2 
4c — 


x2 
2 C Ba a Nan at 
IM | PE =M, A2 GiGoRef angan, 0,0 } : (34) 


This term has quadratic dependence on €, albeit a different one than the total rate (12). 


C. The Bethe-Heitler helicity amplitudes 


In Egs. (31) and (32) we have expressed the leading order HBH rate and interference 


* In this subsection we proceed to 


Ts' 


term in terms of individual BH helicity amplitudes, a 
present the leading order terms of these amplitudes in a readily accessible notation, for 
the special case that the Higgs is at rest in the lab frame. The results below are achieved 
via spinor-helicity techniques; a more comprehensive presentation of the derivation of these 


results is provided in Appendix B. There, explicit results for each spin helicity amplitude 


az, are collected in Eqs. (B2), while the leading order results, obtained by power counting in 
m/mp, are provided in (B3). The compact results below will permit us, in the next section 
of this paper, to study the encoding of the CPV structure in the HBH rate at an analytic 
level. 

Before proceeding, we may first derive a new result concerning the well-studied unpolar- 


2 Using Eqs. (B2) one may show that 


ized BH square amplitude, (Men, |? = Yo... a? 


2 gm E; (k; ` Pi.) — Ey, (ki ` Piz) : _ AG ER. + E? l 
qi (ki + pi_) (Fi Pip) q? (ki: pi )(ki* Piz) 


|MBa, 


(35) 
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This compact expression for the BH rate in the quasi-elastic scattering limit is novel to this 
work. A numerical analysis and validation of the resulting BH differential cross-section is 


provided in Appendix D. 


Now, in the special case that the Higgs is at rest in the lab frame, that is P, = 0, 


Eqs. (B3) for the a; collapse to very simple expressions. Using the leading order results 


(B3) and assuming 0, < 1, for each branch one finds 


an = —lasp)" 


0 2/2 Ya Y— 1 1 
TP Lage. IFE 
wap VART 


d Yeti 


y 


ple t gbe t (36) 
1+7702 1+2 J’ 


1 
21 12 


ap = —lan) S0, 


Here 0+ and H4 are respectively the polar and azimuthal angles defined with respect to the 


branch parent photon momentum, as shown in Fig. 6, and yz = E4/m > 1. We also assume 


that the O(@) terms do not cancel completely. The latter may occur on the phase space 
slice 740, = y_0_ and | — ¢,| = 7, corresponding to minimal |q?|. Excellent numerical 
validation of this expansion implies that the neighborhood of this phase space slice, on which 
the expansion fails, is actually of negligibly small measure. Finally, we may approximate q? 


by 


= @ = m? (7203 +720 + 24 749-04 cos(ġ - 44) | 


Eqs. (36) and (37), when combined with Eq. (28), form one of the central results of 
this paper: a compact form of the HBH rate, which is both numerically stable, and whose 
structure may now be studied analytically. For example, it is now manifest that HBH polar 


angular structure is dominantly controlled just by the k - p,/m? ~ 1+ 404 denominators, 


which produce a peak near 6474 ~ 1. More importantly, we see that the aj, 2, amplitudes, 
that control the helicity interference terms, contain phases which are the leptonic azimuthal 
orientations, ¢;,.. This non-trivial result, when combined with Eq. (32) shows us that, at 
leading order, € manifests as a phase shift in the relative azimuthal orientations of leptons 


with different parent photons. 
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FIG. 6. Definitions of local spherical polar angles {6;,.,¢;,.}. Note in particular that azimuthal 


angles, Hi, are positively oriented with respect to their parent photons, and are defined with respect 


to an azimuth & common to both branches. Polar angles are defined with respect to parent photon 


momentum (black dotted). 


V. SENSITIVITY TO CPV 


In this section we assess the potential sensitivity to CP violation. To do this, we consider 
a number of CPV sensitive observables and propose several sets of kinematic cuts that can 
enhance the CPV signal. These sets of cuts require that the lepton-lepton opening angle 


(13), A, can be resolved, as well as the two photon-lepton angles, 6+. In the following we 


mostly consider angular resolution cuts of the form 


dee j Ox > dent : (38) 


Following Sec. III B, we apply an angular resolution cut Ocu = 1074, which is at the edge of 
what may be possible with present detectors, and a looser, futuristic Ocu = 1075, intended 
to show that very large CPV effects are possible in principle. 

We emphasize that as our measurement is a novel, challenging one, our goal here is not 
to conduct a full collider analysis including backgrounds. Rather our aim is to identify the 
types of observables that can probe CP violation in h > yy, and estimate how well they do 
under ideal circumstances: a high efficiency in reconstructing conversions and a signal rich 
channel. This is in the anticipation that such circumstances might materialize in a future 


LHC running or at a Higgs factory. 
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A. Differential scattering rate 


The HBH differential scattering rate for the full 3 — 6 process (Higgs plus two nuclei to 


two nuclei and two e*e~ pairs) in the lab frame is 
dI x IM |?dITp,_.+dU gu, dgu, 


|M]? dpi, ËP! dk; 
x EME Menlo a e A — — fa) , 
Yy Li=1,2 ta i=1,2 i 


Q= 


(39) 
with |M]? given by Eq. (28). Integrating over the out-going nuclear momenta and all other 


momenta in delta functions, we obtain in the limit |q| « M 


Mp, 


dee 
~ Qn)? 4M4E2 


MP] TI pelan, Jamu ar, 0, , (40) 
i=1,2 


with dQ;, the solid angle for each lepton momentum p;,, and dQ, the solid angle for pho- 
ton ‘1’. The photon labels are extrinsic. We take photon ‘1’, say, as westwards going (if 
Higgs is not at rest one can also take it to be, e.g., the more energetic photon). 

For simplicity, we assume henceforth that the Higgs is at rest in the lab frame, that is 
P, = (m,,0). In the Higgs rest frame the photon angular dependence is isotropic, and 
the integration over dQ, is trivial. Dropping the prefactors, the differential scattering rate 
becomes 


i=1,2 


We may now proceed to consider CPV observables. 


B. Global sensitivity to CP violation 


In principle all the information about č £ 0 (or equivalently € £ 0) is encoded in the 
full HBH differential distribution. The coefficient ¢ may be determined by a matrix element 
method [50-55], as long as backgrounds can be kept under control. Estimating the full 
power of the matrix element method is beyond the scope of this work. 


To test the sensitivity of HBH to č we instead introduce a parameter 
_ J \d'(sm)/dps — dI (c, é)/dps| dps 
E T(SM) | 


ZB (c, Č) (42) 
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Pu > 10°, Co = 0.09 — 
Pu > 10-4, Co = 0.06 


E/T 


FIG. 7. Top panel: Zp in the c— č plane, with the SM point at (c, č) = (—0.81, 0). Bottom panel: 


Z% as a function of € = tan~!(é/c). The scatter of the data points is a numerical artifact. 


The parameter Zg can be thought of as a proxy for the sensitivity of the matrix element 
method, once one integrates over the full phase space. In the top panel of Fig. 7 we show 
the value of Zg in the (c,€) plane. There we see that the deviation from the SM is mostly 
due to the c?+ @ enhancement of the h — yy rate, which need not arise from CP violation. 
Such an enhancement is best detected by measuring the overall h — yy rate, and not using 


our method. 


To assess the sensitivity to CP violation alone, we should restrict our attention to the 
circular contour c? + €? = cm, on which the total HBH rate matches the SM rate for any £. 


This contour is shown as a white circle in the top panel of Fig. 7. To this end, we define a 
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second quantity 


f |at(0)/ars — dÎ (£) /dps| dps 
E aa (0) 


Z3(f) = ZB(8) , (43) 


where dI'(€) is the differential HBH rate with the enhancement of the total h > yy rate 
factored out, 


Tegas- CAP) gece) aps. (44) 
CSM 


Note that in general Zg € (0,00), but Z$ € (0, 2]. 

The bottom panel in Fig. 7 shows the value of Z% as a function of €. The sinusoidal 
dependence of Z% on € is not unexpected. For instance, in the A — yy toy example we 
considered in Sec. IC, in which the angle between linear polarizations is measured, one has 


from Eq. (12) 


ae ca 4 
ZAO) == | [o0s*(o-+ €) — colo) = “Ising | (45) 


h>yy 


By comparing this toy system with HBH, it is therefore natural to deduce that Z$ (£)/| sin €| 
provides a measure of the average CPV signal size obtainable via the matrix element method, 
for any c,¢é. That is, this measure is independent of € and the overall normalization c? + @. 
The fit in Fig. 7 suggests this CPV signal is O(10%) for an angular resolution cut ee > 107°, 
and O(5%) for Ore > 1074. Hence, for the most pessimistic case that there is no deviation 
from the SM in the total h > yy rate the detection of NP from the full matrix element 


method will be challenging, even if there is large CPV component in h > yy. 


C. Differential azimuthal scattering rate 


Let us also consider the sensitivity of an experiment in which just one relative azimuthal 
angle — the difference of the azimuthal angles between two opposite branch leptons — is 
reconstructed, such as the experiment described in Fig. 1. From Eqs. (28) and (36) we saw 
that € manifests as a phase shift in the relative azimuthal orientation of leptons on different 
branches. It is instructive to write down this manifestation explicitly. Let us transform from 


the azimuthal coordinates ¢;, to the coordinates 


P = hi, + Px 3 Ei = Qi — Vins (46) 
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and choose 61, = 0, without loss of generality. With this choice y and £1 € [0, 27), while 
E€ E (—27, 27). For the case that the Higgs is at rest, we find from Eqs. (31) and (36) that 
the HBH square amplitude, and thus the leading order differential scattering rate, takes the 


generic form 
1 
dV /dps = ad e + 2 bj cose; + > Cp COS (NpE1 + MeE2 + 2Y + 2E) | , (47) 
112 j k 


where a, bj and cy, are real functions of yi, and 6;, — they span the energy, polar angle 


phase space, denoted hereafter by PS,» — but are independent of the azimuthal structure, 
and nz, My are positive integers that satisfy ng + mx = 0,1, or 2. From Eq. (37), one sees 
that the 1/q* factors depend on both the azimuthal and polar angles, but in a way such 
that qf = n(1 + (cose)?, where ¢ < 1. Hence 1/q* may be expanded in a power series of 
(cose < 1. Integrating over the £ acoplanarity angles, one may then show that, with respect 
to the azimuthal structure, only constant or cos(2y + 2) terms survive. That is, the leading 
order differential scattering rate 


dV 


— (212 
ie + 2)| Aye + By cos(2p + 26) (48) 


in which the y,@ subscript denotes exclusive dependence on the energy polar angle phase 
space, PS, 9. 

The results (47) or (48) show that £-dependence and y, 0-dependence factorize. Hence € 
appears only in the azimuthal structure as a phase shift. Specifically, the CPV parameter £ 


manifests itself in the inter-branch azimuthal structure in the differential rate 


dv 


ro = (7 + č) Aaea + Bro) ps, , cos(2p + 2£)| , (49) 


which we have now shown is oscillatory with respect to y at leading order. Note that we 
could have chosen y to be any of the four inter-branch angles y1,2, = di, + ¢2,, where 
a, B = +, of which three are linearly independent. Results similar to Eq. (48) follow with 
appropriate replacements. 

In the right panel of Fig. 8 we show the differential distribution dI /dọy for HBH events 
(including a loose cut on the polar angles, at 10~°) for two values of £. It is evident that the 
oscillation amplitude, (B) /(A), is small — approximately ~ 2% — when averaged over all of 


the phase space PS, 9. 
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The key question we now wish to address is whether such a small oscillation amplitude 
is because of small, O(1%) oscillations, or whether there is an O(1%) part of the phase 
space where deviations from the SM are O(1). In the language of Eq. (48), this question 
can be rephrased in a precise manner: Is (B)y/(A)u small for all U c PS, or does there 
exist U C PS such that (B)y/(A)u is O(1). The latter possibility is phenomenologically 
preferred, since it permits the extraction of an O(1) CPV signal on U, which would scale 
better with increasing statistics. 

To address this question, let us begin by examining the coplanar limit. In this limit the 
acoplanarity angles e, » are both zero, and we have from Eq. (47) 


dv 


dy dPS- g = Ayo + Br cos(2& + 2) . (50) 
Y 


In the coplanar limit the size of the modulation is given by 


BN — Ri Ng Vig bip Ti OP (51) 
H9 i=1,2 (1 T (ALAN Ot pa 20?) Ar Rae On + Yiz hi)? | 


where R; = 294, y./ (92, #92 )- The ratio B/A" is small when 70 ~ 1, i.e. near the peak 
of the square matrix element, but B°°/A°° — 1, for y0 > 1 and y+ not much bigger than 


y+. An example is shown in Fig. 8, where E,, and 6;, are held fixed such that 6;, > m/E:, 


and y+ ~ y_. In this slice of phase space the azimuthal oscillation amplitude is large and 
there is a strong sensitivity to CPV. This shows that regions of phase space with large CPV 


signals exist. 


D. CPV enhancing cuts 


We now use the results from Sec. IV B to design kinematic cuts that enhance the sensitiv- 
ity of dI /dọ to CP violation. That is, we seek to enhance the ratio (6) /(A) for a particular 
subset of the HBH event sample. The cuts we propose fall into two classes: those which 
are placed on the kinematics of the whole event; and those which are placed independently 
on individual photons and their daughter leptons. As they use all of the information in 
the event, including correlations among the two photons, one might expect that the former 
produce better CPV signals compared to the latter. As we shall see, however, both classes 
of cuts perform approximately equally well in enhancing (B) /(A) on their respective phase 


space subregions. 
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FIG. 8. Left: Illustration of O(1) oscillations and phase shifts in the HBH differential rate for 
a sample coplanar kinematic configuration. The azimuthal angle ọ in this slice is defined as in 
Eq. (46). The kinematic configuration is: Ej, = Ej = mp/4, Oi, = 10-4, 6. = 20:, so that 


y+0+ ~ 10 > 1 and y4 = 7_, cf. analysis of Eq. (51). Right: The azimuthal distribution dr /dọy 


for € = 0 and for € = 7/4 with a polar angle cut 6;, > 10”? and ĝe > 1075. The modulation 


amplitude is 2%, but will grow to O(1) once optimization cuts are applied, see Sec. V D. 


1. Cuts on collective kinematics 


The helicity interference terms (32) in the HBH rate arise dominantly from the term 


Relar,,ar, 03,,0,€ 28). We can use this observation to pick only events in which 
Q1,,%1,,2,,%,, IS comparable to the rest of the squared amplitude. With reference to 


Eqs. (33) and (34) we thus introduce several sensitivity parameters Jp 


2 2 
IA lor, 03, +» ana + Aal, (52) 
j#k 
ThASk 
where n = MO, RO, J, with 
Am = AO, 444,02) 9%25, | ? AR = 4Refax,, A15, 02,05, | ? X = Alma, 15, A227] uy 
2 
and Aw RJ = ` Vine kpo + XMR,J : (53) 
i#k 
T#8,P#To 


The first two terms of the denominator in Eq. (52) are simply the € independent parts of 
the HBH squared amplitude (31), while the X,s are various pieces of the interference terms. 
In particular, Tm is the magnitude of the full leading order interference term, while 7r,j 


are respectively the CP-even and CP-odd interference terms. We expect each to be useful 
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gauge of sensitivity to CPV independent of €. For example, cuts on Tr and 7j will enhance 
the azimuthal modulations in the € = 0 and € = 7/4 respectively, and can thus be used to 
distinguish among these. In all cases 7; — 1 (T; — 0) implies full (no) CPV sensitivity. 

For all numerical analysis we use three private Monte Carlo codes that were cross-checked. 
The details on Monte Carlo event generation can be found in Appendix F. Placing a cut 
on Jn produces an event sample with large (B)/(A) ratio. This is shown in Fig. 9 for 
opening angle cut 6 > 10~° and in Fig. 10 for opening angle cut of Oc > 1074. The dr /dọ 
distributions are shown for two choices of CPV parameters, € = 0 (blue histograms) and 
E = 7/4 (red histograms). The fits to the functional form (48) of the dr /dy HBH differential 
rate are shown as solid blue and dot-dashed red lines, respectively. The efficiencies of the 
cuts for the examples shown in Figs. 9 and 10 are ~ 1% for the upper panels and ~ 0.1% 
for the lower panels (the exact values of efficiencies depend on the value of £). The presence 
of € £0 exhibits itself as the expected phase shift in dI /dọ differential rate (48). 

From Figs. 9 and 10 we see that it is possible to select regions of phase space such that 
(B) /(A) ~ O(1). In particular the average modulation (B)/(A) is large (small) for € = 0 
(€ = 7/4) for Tkg and vice versa for Tz, as shown most strikingly in the bottom panels. 
In contrast, the average modulation size does not depend on € for the events selected by 
Tm. This suggests that using several of the parameters Jn simultaneously may optimize the 


sensitivity to CPV further. 


2. Cuts on individual photon conversions 


We now turn to discuss cuts on individual photon branches of the HBH process. These 
cuts are generated by simple sufficiency conditions, that ensure a large CPV signal. 


From Eqs. (31) and (32), in order to ensure that the €-dependent terms are comparable to 
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the full rate, it suffices to require |a}, an| Z las, |? on each branch (i = 1, 2) independently. 


For the case that the Higgs is at rest, using (36), this sufficiency condition is implied by 


S 2 1 for each branch, where 


2 0v0- 
S = 2(1 — cos ô) l m p z (1+4262)(1+ 7262). (54) 


For brevity, we have dropped the Latin branch index, and here 6 = |¢_ — +| — 7 is the 


acoplanarity angle. 
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dI /dọ larbit. units] 
dI /dọ [arbit. units] 
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FIG. 9. The differential HBH rate dr /dọ as a function of the azimuthal angle p between the 
positrons setting the CPV parameter to € = 0 (blue histograms) and € = 7/4 (red histograms). 
The histograms are binned Monte Carlo events with different cuts on parameters m, Eq. (52), 
as denoted in the panels. The solid blue (red dot-dashed) curves are the result of fitting the 
normalized binned events to Eq. (50) for the € = 0 (€ = 7/4) cases with € also floated in the fit. 


The angular resolution cut is dee > Beut = 107°. 


It is notable that this sufficiency condition weights the desirable energy-polar angle regions 
of phase space inversely by the amount of acoplanarity. Specifically, S = 0 in the exact 
coplanar configuration 6 = 0. Of course, S is not a necessary condition for a large CPV 
signal, so relatively coplanar events — i.e. events for which 6 < deut of the event sample 
(dcut = 0.4m or 0.257 for our MC, see App. F) — may in principle significantly contribute to 
large CPV signals in d''/dp. Also note that symmetric conversions, i.e. 7,0, ~ y_6_, are 


more likely to produce a large S, and hence a strong interference effect. 


Fig. 11 shows the results of the Sj 2 > 1 cut for the angular resolution cut (38) beut = 1078 
for € = 0 compared to € = 7/20 (top left) or 7/4 (top right). Note we also constrain 


the acoplanarity 61. € [37/4,57/4], since we do not expect extremely acoplanar events to 
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FIG. 10. The same as in Fig. 9, but with the opening angle cut Ore > becut = 1074. 


encode polarization information. The angular resolution cut alone retains approximately 
40% of the total HBH events on this acoplanar domain; for € = 7/4, the corresponding 
azimuthal distribution is shown in Fig. 8. The addition of the S12 > 1 cut reduces the 
cut efficiency to 2.6%, but unlike Fig. 8, now (B)/(A) ~ O(0.2). The phase shift due to 
non-zero € is clearly visible in both plots, and the value of € extracted from the fits agrees 


with the input values of € = 7/20 or 7/4 respectively. 


Finally, the bottom panels of Fig. 11 displays the acoplanarity on each branch for the 
E = 7/4 HBH events with no S cut compared to the S12 > 1 cuts. We see that the S cut 
mildly favors acoplanar events, as it broadens the acoplanarity distribution and disfavors 
relatively coplanar events. For example, we see that events with acoplanarity 6 < 1% are 
disfavored in the cut distribution, and we also see that the full width at half maximum of 
the acoplanar distribution increases by 50%, from 0.047 to 0.067, under the S > 1 cut. The 
excellent performance of S compared to Tn cuts (see Fig. 12 and Sec. V D3 below) therefore 
suggests that acoplanar events — e.g. with ô & dc — play an important role in encoding 


the CPV signal. 
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FIG. 11. Top panels: The azimuthal distributions dI'/dy for € = 0 (grey histograms) and for 
E = 1/20,7/4 (blue histograms, left and right top panels) with S1» > 1 on the domain 61,2 € 
(37/4,57/4). The solid (dashed) curves denote fits to Eq. (48), with € a free parameter in the 
fits. The bottom left (right) panel shows the acoplanarity distributions, dr /dô¢ġ, 6¢ = (9x — ¢_) 
mod 27 for each photon branch, displayed by blue and gray histograms respectively, with € = 7/4 
and no S cut (Si2 > 1). Note that the scale varies between these two plots. The corresponding 
cumulative acoplanarity distributions cdf(d@) for each branch are denoted by solid black and grey 


dashed lines on each plot. In all panels the angular resolution cut (38) of Oeut = 10”? was applied. 


3. Cut scheme efficiencies 


It remains to determine the efficiency of the above cut schemes. In Fig. 12 we show the 
CPV signal (B) /(A) for the various „n and S schemes as a function of the fraction of the 
total MC event sample, not rejected by combined sensitivity and angular resolution cuts. 
Up to small corrections this fraction is the absolute cut efficiency, see Appendix F for details. 


For comparison, the data point with only angular resolution cuts is also shown for each plot. 
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We see that an increase of the CPV signal by an order of magnitude roughly requires an 
order of magnitude penalty in sample size. Moreover, for high cut efficiencies, the S; scheme 
provides the largest CPV signal. In contrast, in the low efficiency, but higher signal regions, 
the 7,, schemes outperform the S; cuts. To give an idea about the feasibility of this analysis 
at the (HL/HE-)LHC, Table I gives the expected number of events after the application of 
the cuts discussed above. These numbers do not inlcude experimental acceptance efficiencies 
and so will be lower in practice. 

The S; parameter is defined on an individual photon branch, independently of the lepton 
configurations on the other branch, while 7;, is defined on the configuration of whole HBH 
events. We therefore might deduce that the main difference between these two schemes 
is whether they are affected by inter-branch leptonic configurations. If this deduction is 
correct, then the comparative performance of S; versus Jn cuts, shown in Fig. 12, suggest 
that inter-branch configurations become more important for the extraction of larger CPV 


signals, but are not important over most of the CPV sensitive phase space. 


E. Triple products 


Finally, we briefly comment on the possibility to use CP-odd quantities such as triple 
products to directly search for CPV. Since the CP odd terms in |M|? are C even and P 


odd, we could consider two such contractions 


Ty = PIPL, PR P2, (55) 


Vs IL [fbt] o x BR(h > y7) [fb]|Events 
8 20 47 0.24 
14 | 3000 125 94 
33 | 3000 444 333 
100} 3000 1875 1406 


TABLE I. Expected number of events after the application of S or T cuts with Ope > 1074 to 
obtain (B)/(A) ~ 20%. The Higgs production cross section includes the gluon fusion and VBF 


channels only and is taken from [56]. 
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FIG. 12. Comparison of different sensitivity parameter cut schemes, in terms of cut efficiency, 
for various angular resolution cuts. The black data point on each plot denotes the efficiency and 
(B) /(A) for its angular resolution cut alone, with no enhancements from sensitivity parameter 


cuts. 


and 


P 
BS (ePP1 P2- P24 Pi Pp2 po, aa eP! P14 Bey ae Pi Pa ), (56) 


pis | 


where P” = (M,0) is the nucleus momentum, and €??* is shorthand for the Levi-Civita 


contraction €4”?"p,qvT 8c. In terms of scalar triple products, 


T = E1, Pi (Po, X po_) — Bi pi, (po, X po_) +102, 67) 


m = M |p: (pp, Xp») — Pip (pp, X P2) +12]. 


Note, that Tə is not strictly speaking C-even, since it involves the nucleus momentum. 
However, BH conversion on nucleus is the same as on anti-nucleus within our uncertainties. 
This is equivalent to leaving P” unchanged under C transformation. A detailed analysis is 
beyond the scope of our work, but we remark in passing that a straightforward use of the 


S and Jn cuts does not lead to appreciable non-zero value of (7). Further investigation is 
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warranted. 


VI. CONCLUSION 


In this work we have studied how to probe the underlying CP couplings of the Higgs 
to two photons, which undergo Bethe-Heitler conversion on nuclei. We have shown that 
sensitivity to CP violating couplings is possible only if there is interference between con- 
version amplitudes with different photon helicity. Using spinor helicity methods, we have 
computed compact, leading order expressions for these amplitudes, which are novel to this 
work. Our analytical control of the leading order h > yy 2X de full differential scattering 
rate permits us to show that: (i) the differential rate with respect to the relative azimuthal 
angles between leptons with different parent photons is oscillatory; (ii) CPV is encoded as 
a phase shift in such distributions; and (iii) we may construct various sensitivity parameter 
cuts that extract the regions of phase space on which such oscillations — the CPV signal — 
are order unity on average. These analytical results have been confirmed and explored with 
numerical simulations, including a comparison of the relative effectiveness of the different 
sensitivity parameters. 

For simplicity we restricted our numerical and CP sensitivity analysis to the case that 
the Higgs is at rest in the lab frame. This is not the case in the LHC experiments, and this 
assumption needs to be lifted for more realistic studies. We note, in this vein, that Eqs. (B3) 
hold also for the boosted Higgs case, although Eqs. (36) do not. The sensitivity parameter 
cut schemes are expected to work comparatively well in the case that the Higgs is boosted, 
too. On a similar note, we expect these sensitivity cuts to enhance the full matrix element 
method. This method, as characterized by the parameter 7%, appears to be a few times 
more sensitive to CPV than experiments measuring a single relative azimuthal angle. This 
leaves open the possibility of further improved CPV signals, compared to those shown in 
this work. 

We do not expect a h > yy BY de experiment to be straightforward, and further exper- 
imental and theoretical studies are needed in order to see if the methods discussed in this 
paper can be used in practice. From an experimental viewpoint, it needs to be determined 
how well one can reconstruct the electron and positron momenta for opening angles ~ 1074 


or even down to ~ m/m, ~ 105, and whether there are significant rescatterings after BH 
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production. Already from our preliminary studies it is clear that large statistical datasets 
and fine granularity of the detectors will be needed, such as at the proposed HE-LHC, VLHC 
or TLEP [57]. Experimentally, the situation may be more favorable in such a machine, hav- 
ing a larger amount of statistics and better kinematic control of the Higgs. A completely 
independent direction for measuring the CP violating coupling of Higgs to photons — and 
a direction not explored in this paper — would be the use of polarized photon beams at a 
photon collider. From a theoretical viewpoint, it also remains to determine and search for 


other CPV sensitive observables apart from the azimuthal distributions discussed here. 


Lastly, while we focussed on Higgs diphoton decays in this work, most of this analysis 
holds for any pair of correlated photons. They can arise from a resonance, like the Higgs, 
or from scattering events. While the details have to be determined in a case-by-case basis, 


the principles, such as helicity interference, are the same. 
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Appendix A: Spinor-helicity formalism 


In this paper we extensively use the spinor-helicity formalism, in which the sigma ma- 
trices g”, solder null momenta to Weyl spinors, that transform under the spinor irreducible 
representations of the covering group SL(2,C). For a review see e.g. Ref. [58] or [59]. 


For a null momentum k, the associated Weyl spinors Ay are soldered via 


Ko IK K3 kuaka = Arda Ana » (A1) 
in which we have written the dotted and undotted Weyl spinors in the notation 
Gia, k= One, [RY = Orla, Ia, (A2) 


so (kto et) = (Ag) aor" (Aca, (Kor = (Ar) tom (Me). Upper and lower indices are 
related by the usual epsilon contractions, in particular coh = Pn In this notation, the 


sign superscripts denote the helicities of the in-going and out-going states, corresponding to 


kets and bras respectively. In particular, note that |k*)' = (k*| by construction, (k*|o“|¢*) 


do not exist, while (k~|o"|€~) = (€*|o"|kt). For null momenta k and £ we often write the 
bilinears 

(k 0) = (|) = (A) Moda [k 4 = (RTE) = Ona d0)", (A3) 
so (l k) = —(k 4, E k] = —|k 4, (k OU k] = 2k - £ and (k &* = [£ k]. 


In this notation, definite helicity polarization vectors associated to k may have form 


(rž oma 
VER) | 


with r” a reference null momentum. One may show via the Fierz relations that these 


— 


et*(k;r) = + (A4) 


polarization vectors satisfy the axiomatic requirements k-e* = 0, e+ - (e*)* =e*-eF = —1 
and e“ -e= = 0. In the present work, a particularly convenient choice of reference momenta 
is 

T= kz, r2 = ki, (A5) 


where rj» (k12) are reference momenta (photon momenta) corresponding to polarization 


vectors er» defined in (7). We shall always make this choice of reference momenta. 
The dotted and undotted Weyl spinors have explicit representations in terms of the 


momentum components in a particular basis, up to a free choice of phase. For example, for 
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a null momentum k = (k°, kt, k?, k3), a possible choice is 


ki + ik? Tg 
Je k? — k , k <0, 


kt — ik? 
0 3 3 
(vk +k’, 5) 0 Sts 


The components of the upper dotted spinor are immediately specified through the relation 


Or)" = (A6) 


(Ar) = (Ok) Y5 and the lower index spinors through appropriate epsilon contractions. With 
this choice and a choice of basis, all spinor objects — e.g. (k £) — can be evaluated explicitly, 
just as momentum objects may be, e.g. k- 4. Note that for our Weyl spinor phase choice 
(A6) one has (ky k2) = (ko ki] = mp. 

The Dirac spinors of massive particles may also be represented in the spinor-helicity 
formalism, via the application of a light-cone decomposition. For a massive spinor of mo- 
mentum p and mass m, we define an associated null momentum 

2 2 


m 
(PF, or equivalent! K = pp + — 
2p- £ q rP j 2p. £ 


eH (A7) 


where (“ is a null reference momentum. The spinors then decompose similarly as 


up) = Ip") + pale Wp) = P3) ca H 

ap) = (p*| + T W], a?(p) = (57| aah 

up) = I) + tal) v?(p) = |5") - pg 

o(p) = (57| + pg a? (p) = (pt | - cA | (A8) 
One may verify that these spinors satisfy the canonical requirements Tu = 2m, dv = —2m, 


vu = wv = 0 and the completeness relations 9”, wu = p +m and 9”, 40 = p—m. Just 
as for the polarization vectors, one is free to choose the null reference momentum. This 
choice amounts to a choice of ‘gauge’, under which the unpolarized square amplitude must 
be invariant, but the polarized square amplitudes are not. In this work, we shall always 


make the reference momenta choices 


bi = ko and ba, = ky r (A9) 


These are convenient choices for the purposes of extracting the leading order BH helicity 


amplitudes, as below. 
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Finally, together with the several well-known spinor identities, the following identity is 


especially useful for computing terms involving Levi-Civita contractions, 


COPY ging o a = 4i (ame S: ee) . (A10) 


Appendix B: BH spin-helicity amplitudes 


Here we write down explicit expressions for the helicity amplitudes a}. From Eqs. (16) 


and (24) one has (suppressing the branch index) 


+ D a p -ktm E-p, 4m 1, 
0 aa pe 9 gr 
O Ois (EEDU WQEN WFR LK” 
E sal ( k - p- k- p4 ) kop. k- » ABI) 


where u = u(p_) and v = v(p;). Applying the light cone decomposition (A8) with reference 
choices specified in (A4) and (A9), gives the full results for each spin helicity amplitude 


an = (at) 1 |- ZEE Kp) (CQIB) . m [Ep4] (k p+) (C1) 
j 2 PV k- p- [ke] (ep_) k - p4 [k 4] (p_) 
_ miep_] kp») (b_*| OM) , mies (k B+) (p_*| Ql) 
k - p- [k £] (LD4) k + py [k £] (C4) 
Mm KOK Qo) J m (k £) (pH can 
kp (pb) k - p4 (L p+) 
i 1 m (ke) KT MY m? (kpt) (C|) 
meen aa aly eedees tita 
mM Kp ep) UAN Kp (kp— H1 Op") 
k - p- [k 4] [E 5+] (ep_) k - p- [k 4 
i [p+] (k B+) P1 Mp) | (k p+) ad 
k + p4 [k 4] k - py 
az, = Ha) = 1 | [kp_](kb_-) 6-1 Op) | W D+] (k B+) (P| Pp) 
21 Bo v2 k-p- Ia k - p4 [k 4] 
pae (KO (E*| Q|k*) m I D] (k p+) (OI) 
k - p+ CP (£ B+) -p+ [kel Pp (E+) 
mp (AM é D1 Op) 
k - p- [k 4] (€p+) k - p- 
az = (ot) = 1 -72 p-)(k | QIE hm (kp P+) 6 1G) 
22 ji v2 k - p- £p] k - p4 [k 4] 
MURI KPA GTA m (kpt) (C7) elk») 
k - p- [k £] Kp] | k - p4 [¢p_| 
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mep (kp) CK 17 m (kp) Ca LO ur a 
kp QP] kph 1 


As per the main text, we have dropped the photon subscripts, and k, l = ky, ka or ks, kı for 
parent photon 1 and 2 respectively. Squaring these amplitudes, taking traces and summing, 
one obtains the full Bethe-Heitler square amplitude that is obtained by the usual Feynman 
methods. 

We can further extract dominant terms of the BH spin-helicity amplitudes by observing 


that if q; < mp, then (kı k2) >> (ki pi) etc. Moreover, in expressions such as aj, or a3), we 


may discard subdominant O(m?/k; - pjzi) terms. This leads to the following leading order 
results in m?/k; + pjz; and (k; p:)/(ki ko) for the BH spin-helicity amplitudes 


m [ik OPE lkt) (k Oppo lk) 


an = —(Qg9)* & Pool ep kp UP 
PPL Cc (CPM 2 2 cl 28 MNC AR 
12 21 PVK - pa q? V2k pik 4] q? V2k - pa [k 4] 


0 mtk 0 ON) 
PV2k- p_[é PNP) 
(k p_)(ptlo°|k*) E p_](k B_) lo ly) fe [I B+] (k B+) Kp lo By) 


an (an) &- @V2k-p PV2k-p-Ik 4 PV2k pk 4] 
__m2{k Ok E) 
P2k- pil PEB 
aj = -h = 0, ve 


in which we have dropped the photon subscripts, and uniform overall signs or factors of 7; 
k,l = ky, k2 or ko, kı for photon 1 and 2 respectively; and the spinor notation is detailed in 
Appendix A. The parity relations (29) are satisfied as expected. Eqs. (28) and (B3) together 
provide a compact expression of the leading order HBH square amplitude. 

It should be understood that the particular form for the spin helicity amplitudes above 
depends on the choice of reference momenta, because the amplitudes explicitly depend on 
polarization vectors and spinors (see App. A). Moreover, the ability to straightforwardly 
expand the full results to the leading order results depends on a sensible choice of reference 
momenta. In contrast, the full unpolarized BH rate is independent of polarizations and 
spinors, as a result of spinor and polarization vector completeness, and therefore must be 


independent of any such reference momenta choice. 
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Appendix C: Polarization-decomposed HBH rate 


Here we give the explicit results of a Higgs-Bethe-Heitler Feynman type calculation. To 
preempt the loss of numerical precision from large cancellations due to the Ward identity, we 
do not use polarization completeness relations. Rather, we retain the polarization vectors 
explicitly in the HBH rate. In the case that the Higgs is at rest in the lab frame, we simply 
use a Cartesian basis for the polarization vectors, aligning the back-to-back photons with 
the z-axis. That is, in the HBH square amplitude we coherently sum over the polarization 
basis 

e (kı2) = (0,1,0,0) and €%(ki2) = (0,0, 1,0) . (C1) 
The squared matrix element for the process h + Ni + No > y(ki)y(ke) + Ni + No > 
4+ Ni + N} is given by 
IM)? = Ar APS let Eb (ko) (ku eG (Hea) ey (kan eG (ken) bo (ka) ef" (ko) | BH] BH3 
Kg 
(C2) 


where the tensor A“”” is 
ae =C (kı ` ko go = kok.) + če ky hog . (C3) 


The BH, factors are the polarized Bethe-Heitler squared amplitudes for a photon 2, including 


form factor contributions, these are in general 


8G 
BH” = E Ak -p Nk) (QE_E, —m?—p,-p_) B® 
EOF || ( p 1 p+) ( t P+°P ) k,p-— 
8E, 
T “ape pa) (k PB, + (k - p4) (4E? + ae + (pis p-)| 
PeP? 
“Skep Pkn E 


(k-p_)(k- p+) (- + 2(E_(k- p_) + Ex(k- p+))] BE, 


+ 2(k-p_)(k- py)(¢? — 4E_E, Be" 


P+,P— 


k- p-)(k- p+) [(k- p- +k- ps)? + PE g” 


EI 


— 4(k - p_)(k - p4) (2E Beat pe) eh, (C4) 


where Ff? = (“kb + (ks — 2690 and Be, = (145 + Gs. The photon momentum is denoted 


by k and the lepton momenta by pi, p_. M and m are the masses of the nucleus and lepton 


Al 


respectively. Assuming that the Higgs is at rest in the lab frame and that the photon is in 
the z-direction, as in Eq. (C1), the previous expression simplifies to 
a 26G {= eg San TN j (Gu us ptt) (Fue i wE) 


q (k - p_)(k- p+) k:p- k:p k-p k:p 
(C5) 


where, in the second term, we expanded terms of the form ,/4E? + q° to leading order 


in g'/E?. We have checked that the helicity formalism results and the Feynman diagram 
calculation results for HBH rate agree. 
Finally, the unpolarized BH rate can be obtained from Eq. (C5) simply by averaging over 


the photon polarization as follows 
(MP) = 1/2 (BH" + BH”) . (C6) 


It is instructive to use the polarization vector completeness relation to obtain an expression 
in terms of Lorentz dot products. In this case, starting with Eq. (C4), the BH rate is given 
by 


(E2 + E3) , m? (k-p_E_ -k:p E) 


(MP) -egf -2 


q? k- p- k- p4 q (k-p_)? (k - p4)? 
me (k:p-+k:p)  k-p_k-py _ (k:p-)?+ (k:p) 
P (k:p) (hope) (k-p-) (k- p4) q (k - p—) (k - p+) 
k- p- +k: p4 
m C7 
nelin) (an 


where the leading terms — the terms on the first line — reproduce Eq. (35). 


Appendix D: Numerical simulations of Bethe-Heitler conversion 


In this appendix we present numerical evaluations of several differential BH rates. The 
numerics were done in two ways: by numerically integrating the full tree-level analytical 
results — i.e. the BH rate arising from Eqs. (B2) or (C4) with appropriate integration 
measures — using the CUBA library [60]; and with a Monte Carlo (MC) code developed 
privately (the details are given in Appendix F). Fig. 13 shows the differential distribution 
of the positron energy fraction E} = E,/E,. For efficiency the MC simulation (blue binned 


histogram) is generated with a cut on the difference of electron and positron azimuthal 
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FIG. 13. Spectrum of the positron energy E+ = E/E}. No opening angle cut was applied in 
the left hand figure and an opening angle cut of 1074 was applied in the right hand one. The his- 
tograms were created with MC events and the solid curves are results of numerically integrating the 
differential cross section. The dashed curve in both figures is the result of numerically integrating 


the differential cross section over the entire range of 6¢ as opposed over the range (0.67, 1.47]. 


angles, d6 = (6, —¢_) mod 27 € (0.67, 1.47] (see App. F). This agrees with the numerical 


integration (red line) for the same cuts on 06. 


Fig. 13 (right) shows the positron energy distribution after applying the opening angle 
cut of A > 10”! on the angle between et and e~ momenta. The MC agrees with the 
full numerical integration of the BH rate even though the 06 cut is still applied in the 
generation of the events . Fig. 13 demonstrates that the asymmetric configurations, where 
one of the two leptons carries the larger part of the photon energies, are the more probable 


ones, especially for non-zero opening angles. 


Fig. 14 shows the positron polar angle distribution and demonstrates the combined effect 
of the Si nuclear form factor and the smaller available phase space that suppress very small 
momentum transfers and thus very small polar angles. The peak is at ~ m/E ~ 10~° both 
for the distribution without a cut on the ete~ opening angle 0x, Fig. 14 (left), and for the 
case where 6%, > 1074 is imposed, Fig. 14 (right). This cut also results in an additional 
peak in the distribution, cf. Fig. 14 (right). 

In Fig. 15, the distribution in the relative azimuthal angle 0 of the two leptons is shown. 
The majority of the events are close to the coplanar configuration, where the photon and the 
two lepton momenta all lie in the same plane. However, it is noteworthy that approximately 


40% of BH events have acoplanarity of ~ 5% or more. 
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FIG. 14. Polar angle distribution of the leptons. No opening angle cut was applied in the left hand 
figure and an opening angle cut of 1074 was applied in the right hand one. The histograms were 
created with MC events and the solid curves are results of numerically integrating the differential 
rate expression. The small bump in the right hand figure (~ 1074) is a result of applying an 


opening angle cut. Its location is a function of the cut. 


Norm. events/bin 


P(ôġ/r € |1- A,1+4A)]) 


0 
0.75 0.8 0.85 0.9 0.95 1 


FIG. 15. Left: Distribution of the azimuthal angle between the leptons. The histograms were 
created with MC events and the solid curves are results of numerically integrating the differential 
rate expression. Right: The cumulative distribution function of the relative azimuthal orientation, 


P(é6¢/n € [1 — A, 1 + AJ), from numerical integration. 


Appendix E: Analysis for qq > yy 


Performing the measurement proposed in this work faces two main challenges: first, 
resolving and reconstructing the electron and positron directions; and second identifying a 
background-poor sample of events with Higgs decaying to diphotons. Regarding the first 


challenge, one might simply ask how well and with what efficiency can the LHC or a future 
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collider detector reconstruct the details of photon conversion. To do so, the experimental 
collaborations may wish to test the polarization structure of a standard model (non-Higgs) 
amplitude. To demonstrate that there is a non-trivial structure to be measured in SM 
conversions, we briefly analyze here the leading production of diphotons at the LHC. 
The dominant diphoton production (and dominant background for Higgs to photons 
events) is gg — yy scattering. This has tree-level spin-helicity amplitude 
Posh kı, Ma 
(Mpa| 422 = + u channel , (E1) 
Pr k2, Ag 
in which r = 1,2 (s = 1,2) is the spin of the (anti)-fermion, and A12 = + are the usual 
photon helicities. The two photons then convert in the tracker. The background rate can 
be written in a form similar to Eq. (28). That is, for one fermion species 


IMpal” = 1G) ` | ` [Areh na a22, ? . (E2) 


T,S T1,81,72,82 A1,A2 


As is well-known, the only non-zero independent amplitudes are (Mpc|i» and [Mga]. 
One finds with our usual choice of reference momenta, and the light cone decomposition 


(A8) for the quark momenta P4, 


A,(ko,k1)  A* (kı, ko) 


[Mae]iz = (Mescla) = Pek Peka ’ 
[Msa] = (lMse|x )” = Ahi au , (E3) 
ay ge E 
Alk, 0) = Q3 P. ael - k] | 
with Qy the fermion electric charge. The square amplitude simplifies to 
IMacl? = 20192 X |[Msa]i ot, az, + [Mega] a7, a3, , (E4) 


T1,81,72,52 
It is interesting to contrast this with the HBH result. Here the helicity interference arises 
in terms of the form aj,,as,ar, A4, rather than from ay,,a5,,a1,,A5,, as we found for 
HBH. Since af, = (a53)", etc, it follows from the explicit results (36) that this helicity flip 


on branch ‘2’ produces a phase change ¢2, — —@2, in the background interference term, 


compared to the HBH interference terms. 
To compare with the Higgs rest frame HBH rate (48), we assume the quark centre of mass 


frame aligns with the nuclear rest frame. Integrating all over azimuthal structure except 
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wv = Dir, — Ho, (cf. p= bi, + 2, in Eq. (46)), one finds that the background rate has the 


form 


ines = ABS + Bee cos(2w) . (E5) 
We see that this differential rate has sinusoidal dependence on the mean azimuthal orienta- 
tion, Y, of the outgoing positrons, with respect to, say, the incoming quarks — the beamline 
— rather than the inter-branch lepton azimuthal orientation, y, as in HBH. Moreover, the 


phase change ¢2, — —@2, ensures background is flat in y. Note also that unlike the HBH 


process, the |Mga]*~ factors ensure this background rate features higher spin waves, so 
that its angular differential structure will differ from the HBH structure, too. 

In summary, the leading-order doubly-converted gg — yy square amplitude is given 
explicitly by Eqs. (E3) and (E4) combined with the BH spin-helicity amplitudes (B3). In 
the gg center of mass frame, the corresponding leading order BH amplitudes are given in 


Eqs. (36). 


Appendix F: Monte Carlo numerical schemes 


To generate Monte Carlo (MC) events for the HBH process, we first generate unpolar- 
ized BH events and then use the von Neumann rejection technique to re-weight the events 
according to the HBH differential rate. For a single HBH event one needs two BH events 
taken from disjoint MC samples. Therefore, we first describe the generation of unpolarized 
BH events. 

The phase space for a single BH event, yN — ete N, is five dimensional. We take 
z-axis to be the incoming photon direction. For conversion of unpolarized photons, the 
kinematics are invariant under overall azimuthal rotations around the z axis. We therefore 
fix the positron azimuthal angle to zero. The remaining four coordinates are chosen to be 
the electron energy fraction E- = E_/E,, two transformed polar angles, tı = log, (#;) and 
to = log,,(0_), and the azimuthal angle of the electron ¢_ (see Fig. 6 for definitions). 

For MC we used two independent private codes. One is written in C and the other in 
C++/Java. To populate the BH phase space we first randomly generate the values for €_, ti » 
and @_ according to either uniform distributions or conveniently chosen initial probability 
density functions (PDFs) and then unweight to obtain the BH event distribution. 

For the C code, the initial PDFs are as follows. For a 60 GeV photon, generating €_ 
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Cc C++/Java 
Parameter Range PDF Range PDF 
E [m/E,,1—m/E,] Uniform |[m/E,,1—m/E,] Uniform 
ty (-7, —2] Uniform (-6, —3| Uniform 
to [—7, —2] Uniform (-6, —3| Uniform 
P— (0.67, 1.47] Lorentzian|  (37/4,57/4| Uniform 


TABLE II. The details on the MC generation of BH events, with phase space variables (1st column) 
for C (C++/Java) generated in the range given in the 2nd (4th) column according to the distribution 


given in the 3rd (5th) column (for details see text). 


according to a uniform distribution results in an efficiency of ~ 70%. We therefore generate 
a uniform distribution of € € [m/E,,1—m/E,]. The transformed polar angle variables, tı 
and ty, are generated according to uniform distributions in a suitable numerical range, see 
Table II. Using the coordinates tı 2 captures the fact that electron and positron distributions 


are sharply peaked around 04 ~ m/E. The BH events are also dominated by kinematic 


configurations that are not too far from the coplanar one. We therefore generate @_ in 
the range (0.67, 1.47], which suffices for our precision. To capture the fact that the BH 
distribution is peaked toward ¢_/a = 1, the C code generates ¢_ according to a Cauchy 
distribution (Lorentzian) with location parameter x, = 7 and scale parameter A = 0.037 [61] 
to improve the efficiency. The C++/Java uses similar initial data, with the exception of 
slightly different ranges and that &— is populated by a uniform distribution on a slightly 
narrowed domain $_ € [37 /4, 52/4] (see Table. II). 


In the next step we unweight the events generated from initial PDFs to obtain the proper 
BH distribution. In the unweighting, the events are rejected with a probability that is 
1—w, with w = (dr ,,/dPS)/ max(dP sn /dPS). The BH MC event sample was validated by 
comparing the generated event distributions to the results of numerical integration of BH 
differential cross-sections as shown in Figs. 13, 14 and 15. Additionally, in Figs. 13 and 14, 
the distributions with a cut on the opening angle between et and e” of 0e > 1074 are shown. 
The MC sample is in excellent agreement with the results of numerical integration. 

In the final step, we convert the generated BH events into MC event samples for the HBH 


process. To do so, two disjoint BH samples were used — one sample per photon branch. The 
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rate for two BH events is given by (dsn, /dPS) (dU gu, /4PS). To obtain the proper HBH even 
rates, we use the standard reweighting technique where events are rejected according to the 


weight 
(dl wen /dPs) kp) 
(dV pn, /dPS) (dpn, /dPS) ? 


where the twist angle between the positrons, wp, is populated by a uniform distribution on 


(0, 27). 


w= 


(F1) 


The two MC codes have been cross tested. In numerics we use 3 x 108 HBH events from 


the C generator and 8 x 10° HBH events from the C++/Java MC generator. 
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